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PREFACE 


The  studv  reported  herein  was  performed  by  personnel  of  the  Ceotet  hnical 
laboratory  (1.1.)  .ind  lnstrument.it  ion  Services  Division  (ISD),  US  Armv  Engineer 
Waterways  Kxperiment  St.it  ion  (WKS),  during  the  period  1  .July  1982  through 
1  lulv  1984.  The  investigation  was  sponsored  by  the  Belvoir  Research  and  De¬ 
velopment  Center  (BRADC),  formerly  Mobility  K<piipment  Research  and  Development 
Command  (MERADCOM) ,  Tort  Belvoir,  Virginia,  under  Project  Order  No.  A2239. 

Tiie  MKRADCOM  Technic.il  Monitor  was  Mr.  A.  H.  C.ranahan. 

The  project  was  conducted  under  the  general  supervision  of  Dr.  W.  F. 
Mareuson  III,  Chief,  Cl,,  and  under  the  direct  supervision  of  Dr.  A.  0. 
Franklin,  Chief,  Earthquake  Engineering  and  Geophysics  Division,  CL.  The  re¬ 
port  was  prepared  bv  Messrs.  R.  F.  Ballard,  Jr.,  and  T.  B.  Kean  II,  Dr.  R.  D. 
lewis,  all  of  Cl.,  and  Mr.  C.  E.  Smith,  Jr.,  ISD. 

COI.  Til  ford  C.  Creel,  CE,  was  Commander  and  Director  of  WES  during  the 
preparation  of  this  report.  Mr.  Fred  R.  Brown  was  Technical  Director. 
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CONVKRSION  FACTORS,  U.  S.  CUSTOMARY  TO  MFTRIC  (SI) 
UNITS  OF  MF-ASURRMENT 


F.  S.  customary  units  of  measurement  used  in  this  report  can  be 
metric  (SI)  units  as  follows: 


Mult  lply _  By 

!«‘t  0.3048 

inches  2.34 

miles  (U.  S.  statute)  1.P09347 

pounds  (force)  4.448222 

pounds  (force)  per  6.894757 

square  inch 


converted  to 

To  Obtain 
metres 
centimetres 
kilometres 
newtons 
ki lopascal s 


3 


ACTIVE  DETECTION  OF  TUNNEI.S  BY  INDUCED  SEISMIC  SPECTRA 

I’AKT  I:  INTRODUCTION 

Bark  gt  mind 

1.  Th»'  piohlem  ol  detect  inn  and  local,  inn  clandestine  subterranean 
.k  t  lvit  it's  has  been  addtessed  for  a  number  ot  years.  Various  geophy s  i  c  a  1 
methods  thought  to  be  applicable  to  the  problem  have  been  evaluated  at  a 
variety  of  sites  and  under  a  variety  of  geologic  conditions  as  described 
by  Ballard  ll^BJ  1.  Many  of  these  techniques  have  been  applied  from  the 
ground  surface  as  well  as  in  boreholes. 

To  date,  only  those  geophysical  techniques  used  on  the  ground 
surface  or  subsurface  geophysical  methods  which  originate  in  boreholes  have 
shown  any  promise  in  the  detection  of  underground  tunnels  or  cavities. 

Tunnel  detection  by  remote  sensing  methods  has  proved  to  be  relatively  ineffec¬ 
tive.  The  use  of  satellite  photography,  infrared  imaging,  etc.,  can  be 
used  to  detect  road  work  or  surface  spoil  areas.  However,  deep-based  tunnel¬ 
ing  activity  has,  thus  far,  eluded  the  state-of-the-art  high-altitude  remote 
sensing  technology.  Such  methods  as  surface  electrical  resistivity,  micro¬ 
gravity,  ground-penetrating  radar,  crosshole  seismic,  and  radar  techniques 
have  proven  to  he  effective  under  the  geological  conditions  which  are  conducive 
to  each  of  the  methods. 

i .  The  thrust  of  this  research  effort  has  been  directed  toward  the 
application  of  crosshole  seismic  techniques.  Recent  work  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  has  led  to  the  development  of 
a  downhole  vibrator  system  which  can  be  accurately  controlled  over  a  frequency 
range  between  SO  to  S00  Hz.  The  device  has  been  successfully  used  in  the 
determination  of  shear-wave  velocities  of  rock  and  soils  between  the  vibrator 
and  receiver.  In  the  course  of  investigating  many  sites  throughout  the 
United  States,  it  has  been  observed  that  earth  materials  exhibit  discrete 
frequency  spectrum  characteristics  which  can  be  attributed  to  many  different 
factors,  some  of  which  would  be  the  granular  structure,  density,  moisture 
content,  state  of  stress,  etc.  This  observation  leads  to  the  suggestion 
that  earth  materials  can  be  characterized  by  their  wave  propagation  spectra 
in  a  manner  analogous  to  optical  spectrographic  analysis  procedures  used 


in  other  fields.  Therefore,  it  an  anomalous  condition  such  as  a  tunnel 
should  appear  between  boreholes  in  the  geophysical  seismic  survey,  its  presence 
should  have  an  effect  on  both  the  arrival  time  of  the  seismic  signal  and 
its  frequency  content.  The  degree  of  change  of  t  fie  total  seismic  signature 
m  frequency  content,  amplitude,  and  arrival  time  is  the  primary  factor 
which  is  addressed  in  this  effort. 


Object  ive 


.  The  bread  objective  of  this  project  is  to  determine  the  modifica¬ 
tions  imposed  on  a  seismic  signal  originating  from  a  repeatable  source  caused 
bv  the  presence  o f  a  tunnel.  It  has  already  been  established  by  Price  (1982) 
and  others  in  previous  studies  that  the  time  of  arrival  of  the  seismic  wave 
train  will  he  lengthened  hv  the  presence  of  an  air-filled  void  (tunnel). 

This  phenomenon  results  from  the  additional  time  the  signal  must  take  to 
travel  around  the  cavity  when  traveling  from  the  source  to  the  receiver. 

These  studies  have  prompted  the  following  specific  objectives:  (a)  identify 
the  seismic  wave  types  that  are  influenced  by  the  presence  of  a  tunnel, 

1:1  determine  whether  or  not  receiver  orientation  (vertical  or  horizontal) 
would  enhance  the  ability  to  detect  a  tunnel,  (c)  to  determine  if  frequency 
i hange  and  amplitude  loss  can  be  directly  associated  with  the  presence  of 
a  tunnel. 


Approach 

.  WKS  has  acquired  two  novel  seismic  sources  which  are  intended 
foi  borehole  use.  An  inhole  vibratory  source  was  developed  primarily  for 
I  fie  determination  of  vertically  polarized  in-situ  shear-wave  velocities, 
kecent  experiments  using  this  device  have  shown  that  it  has  potential  in 
other  applications  as  well.  Accordingly,  t  fie  microprocessor  control  network 
has  been  altered  so  that  the  vibrator  is  now  capable  of  responding  to  program 
mc-d  swept  frequencies,  discrete  frequencies,  and  white  noise  random  exc  ita- 
ti  .  ib  sc  i  va  t  i  "its  using  t  lie  random  excitation  input  mode  in  eon  juct  ion  with 

real-time  fast  Fourier  spectral  analysis,  c ros s - c or  re  1  at i on ,  and  phase  rela¬ 
tions  have  shown  that  definite  spectral  signatures  are  characteristic  of 
different  soil  and  rock  materials.  Additionally,  the  spectral  signature 


car.  In-  drast  icallv  altered  bv  t  hr  presence  of  an  anomalous  condit  ion  such 
as  an  air— tilled  void  between  two  borings  as  theoretically  demonstrated 
bv  i.vtle  and  l*ortn->ft  (1484).  With  these  present  modifications,  even  more 
diagnostic  tunnel  signatures  can  be  obtained.  Typically,  if  tests  are  conduc¬ 
ted  between  two  borings  at  a  number  of  different  elevations  and  the  spectral 
content  recorded,  the  presence  of  a  "seismic  obstruction,"  such  as  a  tunnel 
between  the  buieholes,  will  cause  an  alteration  of  t  fie  frequencies  whose 
wave  lengths  are  compatible  with  that  of  the  tunnel's  geometry.  Not  only 
will  certain  frequencies  he  absorbed  or  shifted  to  other  frequencies  by 
t  lie  presence  of  the  void  (tunnel)  but  the  shear-wave  velocity  will  also 
be  affected  (decreased).  The  second  seismic  source  used  at  the  Idaho  Springs 
test  site  was  an  air  gun  which  was  designed  for  the  generation  of  m i d f reqitenc v 
range  i’-vaves.  In  this  study,  >0  to  r>0()  Hz  were  used.  Since  this  source  is 
c  ssf.it  i.tllv  repeatable ,  as  tests  were  made  in  crosshole  fashion  at  different 
elevations  between  the  boreholes,  alterations  in  the  frequency  spectrum 
could  be  attributed  either  to  material  changes  or  the  presence  of  anomalies 
in  a  manner  similar  to  that  observed  using  the  vibratory  shear-wave  source. 

h .  Tests  conducted  within  the  boreholes  were  designed  so  that  multiple 
receivers  could  be  located  in  more  than  one  boring  at  a  time.  In  the  vicinity 
a  the  tunnel,  tests  were  conducted  with  receivers  at  varying  angles  and  o leva¬ 
ri.  us  ( rom  the  source  so  t ha t  a  two-dimensional  geometry  of  the  tunnel  between 
t  he  !  Telrles  might  he  investigated  by  alterations  in  seismic  signatures. 

St  ope  of  Report 

7.  1  veil  though  numerous  attempts  were  made  to  transmit  a  random  signal 

fr  >r.  the  borehole  vibratory  source  to  receivers  located  in  adjacent  boreholes, 
the  effort  was  unsure esst ul  due  primarily  to  the  relatively  high  concentration 
>!  veit  ic.il  iracturing  which  was  prevalent  at  the  test  site  and  t  lie  low  force 
level  (10  lt.*i  generated  bv  t  fie  vibrator.  As  a  consequence,  this  technique 
was  abandoned  and  emphasis  was  placed  upon  use  of  the  air  gun.  The  primary 
emphasis  of  this  report  will  he  placed  upon  the  presentation  and  discussion 
t  data  obtained  bv  using  t  tie  air  gun  as  a  seismic  source  for  the  crosshole 
tests.  ft  will  also  include  a  brief  description  of  the  vibrator  concept  and 
its  test  because  Its  application  might  prove  successful  at  another  site. 

*  A  table  o'  factors  for  converting  U.  S.  customary  units  of  measurements  to 
metric  (SI)  units  is  given  on  page  d. 
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I’ART  II:  S I TE  DESCRIPTION  AND  IT  SIS  niNDrr.lT.lt 


Idaho  Spr  i  nj^s  Site 


S.  The  test  site  is  located  about  1/4-mile  not t h  of  Idaho  Springs, 
Colorado,  at  an  altitude  ol  approximately  8,000  feet.  Mine  workings  and 
seveial  tunnels  are  owned  and  operated  by  the  Colorado  School  of  Mines  on 
an  expel  imental  basis.  In  the  late  1800';;,  the  workings  were  operated  as 
a  gold  nine  to  exploit  the  Edgar  vein  and  other  less  important  veins.  Rock 
in  the  \ u  mitv  ol  the  tunnels  consist  mainlv  of  intimately  interlayered 
pegmatite,  granite  gneiss,  amphibolite,  biotite  gneiss,  and  microcline  gneiss. 
These  rocks  are  lomplexlv  folded  and  cut  by  pigmatite  dikes  dipping  steeply 
northwest.  A  note  detailed  geologic  description  is  available  t rom  the  C.  S. 
o  .  vi.  T  oiive-.'  ('CCS'  and  the  Colorado  School  of  Mines  (Amuedo  ind  Ivev, 
lw84'.  The  Uelvoii  Research  and  Development  Center  ( BRADC ) ,  formerly  the 
Mob i 1 i t v  Equipment  Research  and  Development  Command  (MERADCOM),  entered 
into  an  agreement  with  the  Colorado  School  of  Mines  to  make  use  of  the  tunnels 
■is  a  test  site.  Several  borings  were  placed  from  a  man-made  plateau  at 
an  elevation  approximately  16U  feet  above  the  centerline  of  one  of  the  tunnels. 
The  borings  were  used  hv  several  contractors  prior  to  the  construction  of 
a  new  adit  which  advanced  the  tunnel  between  several  of  the  borings.  Tests 
were  then  repeated  by  a  number  of  contractors  prior  to  WES  arrival  on  the 
site.  Figure  1  shows  a  surface  layout  of  the  borings  which  were  available 
tiom  the  working  platform.  Most  of  the  borings  were  drilled  to  a  depth 
of  approximately  280  feet  with  a  diameter  of  6  to  7  inches.  The  borings 
were  uncased  except  tor  PVC  pipe  which  existed  in  the  upper  8  to  10  feet 
ot  the  borings.  Tt  had  been  noted  by  one  contractor  that  fracturing  prevented 
the  borings  from  retaining  water,  which  necessitated  remedial  action  by 
grouting  and  redrilling  some  of  the  borings.  This  procedure  was  not  totally 
successful,  so  bentonite  was  used  in  an  attempt  to  seal  the  leaks  and  retain 
more  water.  This  attempt  was  only  partially  successful,  ami  water  could 
only  be  retained  for  a  short  period  of  time  before  drawdown  occurred.  For 
WF.S  purposes,  it  was  necessary  that  water  be  present  only  in  the  air  gun 
seismic  source  borehole  when  conducting  the  crosshole  tests. 
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9.  The  WES  borehole  vibratory  source,  shown  in  Figure  2,  was  developed 
primarily  for  the  production  of  high-quality,  vertically  polarized  shear 
waves.  The  device  is  normally  used  for  the  determination  of  shear-wave 
velocity  through  materials  existing  between  two  or  more  boreholes.  In  prac¬ 
tice,  the  vibrator  is  locked  into  position  within  a  borehole  by  an  inflated 
rubber  bladder.  One  or  more  geophone  receiver  packages  are  locked  into 
position  in  a  similar  manner  in  adjacent  borings  at  the  same  elevation  as 

the  vibrator.  The  borehole  vibrator,  controlled  by  a  power  amplifier  and 
sine  wave  oscillator,  is  then  swept  through  a  frequency  range  of  approxi¬ 
mately  SO  to  SO 0  Hz.  The  receivers  are  constantly  monitored  during  the 

t  requell.  v  -weep,  and  those  t  requeue  ios  which  propagate  well,  as  evidenced  bv 

b.  i  o',  .mip  1  i  t  ode  s ,  are  noted.  A  tone  burst  generator  is  then  used  to  transmit 

a  hurst  of  a  given  number  of  signal  cycles,  usually  eight,  of  one  of  the 
chosen  frequencies.  A  geophone  incorporated  within  the  vibrator  body  is 
used  as  the  reference  standard  to  compare  the  phase  delay  which  is  recorded 
by  the  receivers.  That  difference  in  travel  time,  when  combined  with  the 
distance  source  and  receiver,  is  used  to  determine  the  velocity  of  the  shear- 
wave  within  the  medium.  Using  this  procedure  at  a  number  of  different  test 
sites  throughout  the  continental  United  States,  it  has  been  observed  by 
WES  personnel  that  different  sites  exhibit  markedly  different  preferences 
for  frequencies  which  transmit  well.  In  other  words,  the  earth  materials 
act  as  hand-pass  filters  for  specific  frequencies.  Also,  when  anomalous 
conditions  occur,  such  as  highly  fractured  rock  or  possibly  air-filled  voids, 
a  noticeable  decrease  in  received  high-frequency  signal  content  can  be  ob¬ 
served  is  noted  bv  burr  '  f  1  **H  1 )  . 

10.  Since  the  vibrator  is  an  electromechanical  device  whose  performance 
can  be  controlled  by  electrical  signal  input,  experiments  were  performed 

on  the  WES  reservation  to  determine  the  uniformity  of  signal  output  if  a 
white  noise  random  input  source  were  used  as  the  controlling  factor.  White 
noise  was  chosen  because  it  has  the  characteristic  of  equal  amplitude  content 
in  all  frequencies.  Consequently,  signal  loss  in  any  frequency  band  can 
be  related  to  energy  absorption  within  the  medium.  The  amplitude  spectrum 
of  the  borehole  vibrator  in  the  random  excitation  mode  proved  to  be  essentially 
uniform  over  a  frequency  range  of  80  to  S00  Hz.  This  plot  is  shown  in  Figure  3. 
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Win  If  tins  test  was  being  conducted,  a  receiver  was  placed  at  the  same  depth 
t.’O  led  1  in  an  adjacent  borehole,  approximately  SO  feet  laterally  from 
the  borehole  vibrators  source.  The  amplitude  spectrum  received  at  this 
geophone  is  shown  in  Figure  4.  It  tan  be  seen  from  thts  figure  that  the 
peak  frequency  on  this  curve  is  approximately  22S  Hz.  This  compares  quite 
well  with  other  observations  using  a  swept  sine-wave  frequency;  a  frequency 
ot  221  Hz  exhibited  less  signal  attenuation  or  loss  than  any  other  frequency. 

11.  As  previously  stated,  tests  of  this  type  were  attempted  at  the 
Idaho  Springs,  Colorado  site,  but  the  highly  vertically  fractured  material 
caused  such  severe  signal  losses  that  no  reliable  data  could  be  obtained. 
Several  attempts  were  made  using  lifferent  borings  and  elevations  within 
the  borings  before  it  was  decided  that  further  attempts  should  be  abandoned. 

All  ot  the  crosshole  data  obtained  at  this  site  used  an  underwater  air  gun 
as  the  seismic  source.  These  tests  will  be  described  in  subsequent  paragraphs. 

Air  Cun  Seismic 


12.  The  Bolt  DSOC.I  air  gun  is  t  pneumatically  chargeable,  electrically 
fired  device  (presently  under  study  for  use  as  a  multi-wave  generating  source). 
The  ait  gun  system  consists  of  a  high-pressure  air  or  gas  source  connected 
to  the  gun  through  a  high-pressure  teeder  hose  and  a  battery-powered  firing 
control  box  connected  to  the  gun  through  a  stress  member  reinforced  co-axial 
cable.  The  air  gun  system  is  usually  used  in  offshore  seismic  reflection 
surveys;  therefore,  in  order  to  use  it  as  a  crosshole  survey  source,  several 
modifications  have  been  made.  Because  the  system  is  normally  mounted  on 
a  marine  vessel  equipped  with  high-pressure  compressors,  the  first  modifica¬ 
tion  had  to  be  the  air  pressure  supply  system.  The  system  now  uses  2,500 
ps l  nitrogen  bottles  to  supply  the  required  gas  pressure.  This  decision 
had  two  beneficial  effects;  it  eliminated  (a)  the  need  for  a  large  volume 
compressor  with  its  attendant  logistic  problems,  and  (b!  high  background 
seismic  noise  problems  caused  by  running  of  a  compressor. 

11.  In  addition  to  the  gas  modification,  two  additional  modifications 
were  made  to  the  air  gun;  (a)  a  piezo  crystal  accelerometer  was  added  to 
the  firing  system  by  the  factory  (this  modification  was  needed  to  activate 
the  recording  seismograph),  and  (l>)  a  hydrophone  was  suspended  5  inches 
below  the  gas  exit  orifice  to  sense  a  "zero  time"  for  the  pulse  initiation. 
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1  “  •  The  icu-ivini:  units  used  in  this  experiment  were  two  calibrated 
downhole  tii.ixi.il  geophone  st't  s  (Mark  Products  I.- 1  5 -  fDSWG  >  .  Each  t  ri.ixi.il 
geophone  sot  lonststed  ol  one  vertical  geophone  and  two  horizontal  geophones 
mounted  at  i ight  angles  to  one  another.  These  downhole  geophones  were  handled 
hv  a  stress  memhei  and  connected  to  the  different  recording  devices  through 
an  electrically  shielded  cable.  To  obtain  the  desired  geophone- to-borehole 
sidewall  coupling,  the  downhole  geophones  were  equipped  with  high  pressure 
inf  1  at  able  bladders.  In  order  to  minimize  extraneous  ground  noise  and  vibra¬ 
tions,  these  bladdets  were  connected  through  their  air  lines  to  nitrogen 
bottles  instead  ot  the  usual  small-volume  compressor  generally  used  to  supply 
air  to  t  he  system. 

1  ’> .  In  order  to  build  some  degree  of  redundancy  into  the  program 
and  provide  on-site  "quick  look"  capability,  two  data  collecting  systems 
were  employed.  One  system  was  an  SIE  49R  24-channel  seismograph  mated  to 
a  /-channel  analog  tape  recorder  which  was  used  to  make  analog  recordings 
of  the  wave  trains.  These  recordings  were  later  returned  to  WES  where  they 
were  digitized  and  processed  by  the  Instrumentation  Services  Division. 

(This  procedure  is  described  more  fully  in  the  "Data  Reduction"  section 
of  this  report.)  The  second  system  was  the  EGM.  Geometries  ES1210F  digital 
12-channel  seismograph  which  was  used  to  make  permanent  hard  copies  of  the 
tests  so  that  each  test  would  have  a  uniform  relationship.  All  tests  in 
the  crosshole  series  were  conducted  with  amplifier  gains  pre-set  at  an  optimum 
level  and  unaltered  throughout  the  tests  so  that  all  data  could  be  directly 
compared  on  a  relative  basis.  To  maintain  this  standardization  and  because 
the  SIE  recorder  combination  did  not  have  enhancement  capabilities,  single 
shots  of  the  air  gun  charged  to  500  psi  were  used  with  both  of  the  two  seismic 
recording  systems.  The  sequence  used  to  record  each  event  was  one  shot 
with  the  EG&C  connected  to  the  receivers  and  another  single  shot  in  rapid 
succession  with  SIE  and  tape  recorder  connected  to  receivers. 

lb.  A  variety  of  test  geometries  were  employed  so  that  various  wave 
reactions  to  an  anomalous  condition  could  be  studied.  Tests  were  conducted 
by  repositioning  the  source  and  receivers  at  various  elevations  in  different 
holes.  Figures  5  through  12  show  plan  (a)  and  profile  (b)  views  for  all 
test  conditions  reported  herein.  Tests  were  numbered  consecutively.  Each 
test  number  is  followed  by  a  letter  designating  the  borehole  In  which  the 
receiver  was  fixed.  The  scheme  of  the  data  acquisition  program  intentionally 
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Figure  6.  Test  2  F  schematic 


Figure  8.  Test  1  C  schematic 


(a)  Plan  view 


Figure  10.  Test  4  C  si  tiem.it  ie 


Prot i le  view 


.ill owed  for  (ho  collection  of  more  data  than  would  fie  used  for  analysis. 

In  ('tiler  to  make  this  report  as  concise  as  possible,  only  those  tests  which 
served  to  illustrate  the  methods  used  were  selected  for  presentation. 

17.  The  test  patterns  fall  into  three  basic  types.  The  first  pattern 
is  the  "Common  Depth"  crosshole  test.  This  test  is  used  to  search  first 
lor  changes  in  the  velocity  of  materials  with  depth,  layer  variations  and, 
in  this  case,  the  presence  of  an  anomalous  condition.  In  the  common  depth 
crosshole  test  the  source  and  receiver  are  both  maintained  at  the  same  depth 
during  the  test.  For  this  set  of  tests,  *  he  source  and  receiver  were  raised 
in  1-ft  increments,  beginning  near  the  bottom  of  the  boreholes  at  a  depth 
of  184  feet  and  progressing  upward  to  the  142-foot  depth.  After  this  test, 
the  two  units  were  then  raised  to  130  feet,  120  feet,  and  100  feet,  respec¬ 
tively.  The  following  are  the  common  depth  tests  reported  herein: 

Test  No.  Figure  No. 


IF 

iF 

1C 

AC 


3 

7 

8 
10 


18.  The  second  test  pattern  used  was  a  crosshole  fan  or  tangent  pattern 
where  the  source  was  held  constant  in  each  shot  hole  at  a  depth  of  162  feet. 
(At  the  time  of  testing,  this  was  thought  to  be  the  tunnel's  center.  Later 
information  revealed  the  correct  depth  to  be  about  168  feet.)  The  geophones 
were  advanced  upward  in  3- foot  increments  from  18A  feet,  in  the  same  manner 
as  the  horizontal  test.  The  following  are  the  fan  pattern  tests  reported 
herein: 


Test  No. 
2  F 

2C 


riC 


Figure  No. 
6 
9 
11 


19.  The  third  test  pattern,  designated  Test  13F,  was  the  crossholc 
skew  pattern.  The  source  was  first  held  at  141  feet,  and  the  receiver  lowered 
in  3-foot  increments  until  it  was  9  feet  below  the  source.  From  this  point, 
both  the  source  and  geophone  were  lowered  together  in  3-foot  Increments 
until  the  source  was  at  171  feet  and  the  geophone  was  at  180  feet  (Figure  12). 


2A 


! 

i 

JO.  Tests  conducted  at  Idaho  Springs,  Colorado,  indicated  that  cross- 
hole  seismic  tests  can  provide  extremely  useful  information  which,  when 
properly  interpreted,  can  reveal  the  presence  of  a  tunnel.  The  techniques 
of  data  interpretation  wilt  be  addressed  in  the  subsequent  section. 


PART  III:  TEST  RESULTS 


Data  Reduct  ion 


21.  The  data  which  were  acquired  at  the  Idaho  Springs  test  site  were 
reduced  in  a  variety  of  ways-  The  records  obtained  on  the  EC6G  seismograph 
were  not  recorded  on  tape.  Consequently,  their  use  were  limited  to  optical 
displays.  A  separate  analysis  using  these  records  will  be  discussed  apart 
trom  data  obtained  on  analog  tape. 

22.  That  portion  of  the  field  data  recorded  on  analog  magnetic  tape 
used  an  SI  I  a9-R  seismograph  and  a  RACAI  7  DS  7-channel  FM  tape  recorder.  Tape 
speed  was  1.75  ips  which  provided  a  bandwidth  of  !>.('.  to  1,250  Hz.  The  data 
were  tirst  digitized  at  a  i ate  of  20,000  samples  per  second.  Time  resolution 
.’!  il.OS  nsoi  per  sample  and  a  >  or  respond  ing  Nvquist  frequency  of  10,000  Hz 
were  obtained.  1‘h  l  s  digitizing  rate  far  exceeded  the  frequency  bandwidth  and 
time  resolution  ot  both  the  recording  equipment  and  recorded  signal. 

21.  In  preparation  tor  analysis,  two  types  ot  reduction  were  performed 
on  the  digitized  data.  lime  domain  playbacks  were  produced  and  visually  in¬ 
spected  to  determine  general  signal  characteristics  such  as  time  .if  arrival 
and  amplitude  stand-out.  They  were  also  compared  to  digital  data  taken 
in  the  field  bv  an  EG&G  1210F,  12-channel  seismograph.  After  the  quality 
of  the  data  was  established,  frequency  domain  techniques  were  employed  to 
study  the  propagated  seismic  energy. 

2U .  The  frequency  domain  Btudies  were  performed  using  Fourier  analysis 
techniques.  Fourier  analysis  is  the  analytic  decomposition  of  a  wave  form 
as  a  weighted  series  of  sinusoidal  functions.  This  decomposition  is  done 
via  the  Fast  Fourier  Transform  ( FFT )  algorithm  as  suggested  by  Rabiner  (1975). 
For  the  digitized  record,  consisting  of  N  discrete  values,  the  Fourier 
transform  is  defined  as: 


X  (  k  ) 


X  (  l ) 


(1) 


where 

X(k)  *  the  (complex)  Fourier  transform  coefficient  for  the  k^' 

frequency  component 

X(i)  “  the  i1*1  discrete  value  of  the  input  time  domain  wave  form 
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N  -  t  hr  total  numbet  ot  point  s  in  t  hr  rriord 


(Division  hv  t  hr  total  numhrr  ot  points  in  the  rriord  is  done  to  normalize  the 

»rsu:  t  ito:  trio-.)  The  wave  I  oiri  i.,  t  reh  .  .  de.  onp.iscl  into  .  oil  f  if  lent  s  of 
o  •-  i  ntf.o  1 1  i .  I  i  pel  •.  'lorij.i  1  .  The  «>t  ip  i  n.t  I  time  alien  iv  now  lie  expressed  as  a 
'•nr  . ' !  I  mi  i  i  e  t  •  .  ie  rt  j .  i  on  t  ■•  mi  lit  i :  *  1  it  ‘  (e.  a  .  .  'trip  !  .  ■  ■  .  \  notion  t  i  ,  \  ]  ; 

X<  i  )  =  ^  XU)  exp  (  2  ) 

k  =  1 

Therefore,  the  Fourier  series  coefficients  are  the  amplitudes  of  each 
f requeue v  component  in  the  original  time  series  wave  form.  (Note  that  the 
Fourier  coefficients  are  complex  quantities.  In  this  report  the  modulus  - 
phase  notation  for  complex  numbers  is  used.  The  Fourier  modulus  and  phase  are 
referred  to  as  the  amplitude  and  phase  spectra,  respectively.) 

2 S .  The  Fourier  coefficients  are  also  useful  in  determining  the  amount 
of  energy  a  time  series  has  in  particular  bands  of  frequencies.  From 
Parscval's  theorem  one  has: 


£  I  xm  r  *  £  i x < k )  i 1  n) 

1  =  1  k  =  1 

This  expression  states  that  the  total  squared  amplitude  of  a  time  series 
summed  over  all  time  is  equal  to  the  total  squared  amplitude  of  its  Fourier 
t  r.inst  orni  summed  over  all  frequent  its,  as  shown  hv  Rabiner  (1  ‘l/S).  The  sum¬ 
mation  ot  the  squared  time  series  amplitude  over  all  time  is  referred  to  as 
its  energy.  Consequently,  this  theorem  shows  a  relationship  between  the 
energy  of  a  wave  form  and  the  square  ol  its  amplitude  spectra. 

2b.  If  one  requires  the  energy  content  m  a  band  of  frequencies, 
one  simply  sums  the  squared  amplitude  spectra  of  the  frequencies  in  that 
band.  A  running  summation  of  the  squared  amplitude  spectra  versus  frequency 
will  display  the  total  cumulative  energy  as  a  function  of  frequency.  One 
analysis  conducted  for  this  study  utilized  the  running  summation  method 
of  displaying  the  total  cumulative  energy  (Tests  3F,  SC,  and  13F).  The 
displays  were  normalized  by  dividing  by  the  total  energy.  The  result  is 
fraction  of  total  cumulative  energy  versus  frequency.  This  analysis  is 
referred  to  as  energy  buildup  versus  frequency. 


-  .  i 1  « '  main  I  1 1  .pion.  v  t  n  hn  iquus  usct  in  this  report  wi-if  displuvs 

'  ; '.i  t ^T; p  1  it'll!#'  spin  t  i  i  .in.!  i'iu‘1  v,v  buildup.  Those  two  lei  hniqiics  rrinfnri'c 

<■  >'  '  :!<  r.  Miit  is,  t  hev  ■  mivi'v  similar  intnrm.it  Inn  <1  i  sp !  .i  veil  in  two  dif- 

lotint  w.i-.s.  Mu-  .impliiu.il'  spm  t  r.i  show  t  lit-  il  i  st  r  i  tin  t  inn  of  energy  .is  a 
i'.n.tin  '!  t rcquon.  v ,  Hu'  energy  buildup  displays  tin-  niinu  1  at  i  ve  fraction 
1  t  -t  a  i  energy  as  a  tun  t  i  ■ .  n  of  trequencv.  These  methods  do  stiow  tii.lt  f  re- 
s'a-n.  •  .  'utf.it  an!  i  m  t  gv  dist  t  Unit  ion  varv  as  seismic  energy  is  passed  near 
t  'a  1  nine  '  . 

.M.  A 1 ong  with  the  two  major  techniques  discussed  so  far,  several 
othei  procedures  were  employed  for  analysis  of  the  data.  They  included 
c r oss - c or r e 1  at i on ,  phase  spectrum  analysis,  amplitude  envelope  studies, 
and  energy  buildup  versus  time.  h.u  li  ot  these  techniques  hud  its  merits 
and  drawbac  ks . 

29.  Cross-correlation  is  a  method  of  determining  the  degree  of  similar¬ 
ity  between  two  time  series  (x  and  y)  when  a  time  lag  is  applied  to  one. 

The  c ros s - cor  re  1  at l on  is  defined  by  Claerbout  (  1976  )  as: 

N 

s(i)”~  TL  x  ( k  )  y  (  k  +  i  )  (O 

k  *  l 

whe  re 

s(i)  =  the  cross-correlation  at  a  time  lag  of  i  units. 

Note  that  if  the  y( )  time  series  were  a  t  iroe  delayed  version  ol  the  x()  time 
series,  then  t tie  cross-correlation,  s(),  would  have  a  maximum  when  the  time 
lag  was  equal  to  the  delay  time  between  v()  and  x().  Hence,  the  cross- 
correlation  contains,  among  other  things,  the  time  delay  information.  The 
Fourier  transform  of  the  cross-correlation  is  known  as  the  c ross- spec t rum 
and  is  the  amplitude  of  common  frequency  components  between  two  signals. 

Both  of  these  operations  were  performed  on  the  data  from  the  horizontal 
geophones;  however,  neither  produced  any  usable  insight  into  the  tunnel 
identification  problem. 

30.  Another  scheme  of  tunnel  identification  was  to  explore  the  possibil 
ity  that  the  tunnel  produced  a  dispersive  effect.  That  is,  the  group  velocity 
the  velocity  at  which  the  seismic  energy  travels,  would  be  different  for 
different  frequency  components.  This  dispersive  effect  might  be  detected 

in  two  ways;  first,  as  a  phaae  shift  in  particular  frequency  ranges  or  a 
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i  h.inge  in  .iriiv.il  time  u|  tin-  m.iin  amplitude  envelope.  To  study  this  effect, 
the  phase  spin  t  i  .1  and  envelope  amplitudes  were  used.  However,  the  phase  spei 
tin":  wa  ••  iot  ,  .  ui  i  ,ti- 1  i  t  l>e  part  t<  nf.tr  lv  iim'I  ill  in  the  examined  f  requeue  jes 
am!  the  ,•  i  veil  ••  i  na  1  -  t  .  ■  n.  >  i  :.r  rat  i  •  > .  1  he  ,  hnaeteristies  <>t  the  amp  1  i  tilde 

envelope  may  he  explored  in  the  future  when  a  higher  energy  seismic  source 
is  used.  It  is  possible  that  the  distances  between  boreholes  may  not  have 
been  enough  tor  dispeision  to  he  detected.  Also,  any  noise  present  in  the 
t  ime  series  produces  even  more  noise  in  the  phase  spectrum.  Hence,  the 
dispeisive  effects  mav  have  been  hidden  within  this  procedure. 

11.  Another  method  to  detect  deviations  in  group  ve 1  itv  near  the 
tunnel  is  the  energy  buildup  verso-,  time  ui.ilv.is.  Here  the  traction  of 
total  energv  versus  time  delay  is  plotted.  Although  this  technique  shows 
some  information,  it  was  not  exploited.  One  of  the  problems  with  this  method 
is  noise.  In  this  case,  t  tie  noise  adds  to  the  energy  buildup  even  though 
the  seismic  energy  arrival  does  stand  out  over  the  background. 

D  a  t  a  A  n  a  1  y  s  i  s 


EG4»< .  records 

(2.  The  data  analysis  will  be  addressed  in  the  following  sequence: 
First,  the  time  domain  signals  recorded  on  the  EG6G  seismograph  will  be 
discussed,  followed  bv  the  d i g i t a  1 1 y- processed  data  obtained  on  the  SIE 
analog  tape  recorder.  The  digital  analysis  was  guided  by  the  original  time 
domain  records  obtained  on  the  EG&C  unit.  Even  though  all  data  were  digitally 
processed,  it  was  felt  that  onlv  those  tests  which  served  to  illustrate 
the  analysis  techniques  would  be  used.  Gonsequen t l y ,  Tests  3F,  3C ,  and 
1  IF  wer  subjected  to  frequency  domain  analysis.  Test  13F  was  also  subjected 
to  further  time  domain  analysis  and  will  also  be  included. 

33.  Analysis  was  performed  on  the  EG&G  data  by  first  re-copying  and 
arranging  the  records  in  order  of  increasing  depth.  Then,  arrival  times, 
travel  distances,  and  wave  velocities  were  computed  and  tabulated. 

34.  Tests  IF,  2 F,  and  3F  were  conducted  at  two  different  hole  spacings. 
Tests  IF  and  2F  were  between  holes  G  and  F  over  a  distance  of  about  33  feet. 
Test  3F  was  conducted  between  holes  H  and  F  over  a  distance  of  63  feet. 

35.  Test  IF.  Test  IF  was  a  cornnon  depth  crosshole  test  illustrated 
In  Figure  5.  The  time  domain  signals  recorded  at  each  depth  are  shown  in 
Figure  13.  Table  1  is  a  tabulation  of  P-wave  arrival  times  and  associated 
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v»'  1  oi  1 t  n-  s  .  Test  IK  shown!  imisislfnt  atriv.il  t  lints  of  3 .  3  msec  until  a 
depth  o!  Inn  (fit  was  tcaihed.  At  this  time,  the  presence  of  the  tunnel 
i an  he  seen.  Within  the  tunnel  zone,  t  ime  delays  of  about  1  msec  or  a  total 
lapsed  time  of  4 . 3  msec  were  prevalent.  As  test  inn  progressed  below  the 
bottom  of  the  t  nine!,  t  lie  I'-w.ivc  arrival,  omt  avail),  occurred  at  3 .  3  msec 
until  a  more  competent  zone  was  reached  at  a  depth  of  184  feet,  at  which 
time  tlie  I’-wave  arrival  decreased  to  .’.3  msec.  Velocities  are  plotted  as 
a  t unc t  ion  of  depth  in  Figure  14.  Thev  averaged  about  9,400  fps  above  the 
tunnel  to  a  low  of  7,  300  fps  when  the  tunnel  was  encountered.  The  competent 
zone,  at  a  depth  of  184  feet,  exhibited  a  P-wave  velocity  of  11,000  to  13,000 
tps.  Observing  the  time  domain  signal,  one  can  see  an  apparent  phase  shift 
and  the  signal  is  severely  attenuated  because  of  the  presence  of  the  tunnel. 

ih.  Te_st  .  The  plan  and  profile  views  of  Test  2F  are  shown  in 
Figure  h.  The  test  was  conducted  by  holding  the  source  constant  in  boring 
C.  at  a  depth  of  lb.’  feet  and  moving  the  receiver  in  3-foot  increments  in 
hole  F.  This  resulted  in  a  fan-type  pattern.  Observing  the  time  domain 
plot  in  Figure  13,  the  first  arrivals  of  the  P-wavc  were  determined  and  tabu¬ 
lated  in  fable  .’ .  .’be  shortest  distance  between  source  and  receiver  would, 

<•{  course,  he  it  a  depth  >P  lb.’  feet.  As  the  receiver  moved  in  1-foot  incre¬ 
ments  above  and  below  this  point,  the  distanees  to  the  receiver  increased. 
I'sing  these  distances  divided  by  the  appropriate  arrival  times,  apparent 
velocities  were  computed  and  plotted  as  a  function  of  receiver  depth  in 
Figure  It).  This  presentation  yielded  a  rather  explicit  "picture"  of  the 
tunnel  because  of  the  radical  decrease  in  the  velocity  when  the  wave  travel 
path  was  interrupted.  The  test  conducted  at  a  depth  of  1 fo 2  and  163  feet 
showed  marked  decreases  in  signal  amplitude.  However,  an  increase  in  ampli¬ 
tude  was  noted  from  a  depth  of  167  feet  to  184  feet  .  This  apparent  increase 
in  amplitude  could  not  be  explained. 

37.  Test  3 F .  Test  3F  was  a  common  depth  crosshole  test  with  the 
seismic  source  placed  in  hole  H  and  receivers  simultaneously  in  holes  F 
and  E.  Plan  and  profile  views  are  shown  m  Figure  7.  The  data  acquired 
in  hole  F.  were  not  of  sufficient  quality  to  warrant  m-depth  analysis. 

This  was  because  the  decision  was  made  to  keep  all  gain  controls  at  the 
same  setting  so  that  relative  comparisons  between  different  locations  could 
be  made  at  any  time.  As  a  consequence,  the  amplitudes  were  too  low  for 
accurate  measurement  of  first  arrival  times.  Figure  17  is  the  EC6C  time 
domain  signal  recorded  in  hole  F.  The  tabulated  velocities  and  first 


DEPTH,  FT 


TEST  1  *  F 


0 

100  j-  - 

i 

I 

r 

t 

i 

120  (- 


140  \- 


160  ^ 


180 


280  V 


220  |- 


4000 


VELOCITY.  FP8 
8000  1 2000 


16000 


"T - V  T — 


<3> 


«> 


<}> 

❖ 

❖ 

$ 

«> 


$ 


<$> 


a 


20000 
- 1 


240 


FiK'ire  1  ■’* .  Velocity  viTH'is  Depth,  Test  IK 


1.’ 


t  *  *;  t 


TANGENT  TEST  OR  FAN  OUT 


1 


Ml  i  f  IVf  M  ’>  ’ 

d#pth,  rt 


TEST  2  *  F 


COMMON  Dl PTH  CROSS  HOLE  TEST 


»{  r(.ivCw 


d«ptK,  H 


“  V  4  •  .  *-♦  ^  >  V  >  '  *  * 


A  A  ,A. 

;  \  v 


.  .*9+.,  ■  mH.  iH 


.  !>*  i 


£4  A 


llli  A*  1***^  ,  _ 

■  ,  !  ~  *  v<  r'^  ^  *- 


.  A 

1/  \. 


'  •  t . •  •*■  ; v '  ■ •  •  *, • 


Jk%  ‘A  tm.  m 

■  v.  V  >' 


—  W"*“V 


^  JlV  /+*■ .  ^ 


A  /A  iiw 

1  . 1  \  * 


>v  ^ 

■  •  Vv  > 


-»  1 


!.,M  /A  ,  ,  A-,  A  <^V  A  A,  «k  4»  . 

v<  •  ■  1  v'  v  v  V  ‘  W  *’  "  *”  ’  v  '  “* 


V  y  \t"r  u  v**  '*• 


li»,.4,i..,|CUf 

< 


U  fl  *  » 

«i  i.i  i vi  n  *c’ 


COMMON  DEPTH  CROSS  HOLE  TEST 


d*ptf\  H 

ICO 


i  6  > 


1 


'ft* 


1  7  ? 


’  75 


U8 


i«  r 


Tfl4 


v7 


k  -*i  A  it  A  A  A  a  *, 

v  ^  V  v 


’  V'  •»  * 
1  v’ 


«  A 

*-4  rf**  H+  *  4  *  44  < 


‘  r  •"V 


'-.I 


ft*" 

M'*t-41 

v,,  ^  V,  t*'  -u+l  *  •  * '  *  n; 

‘ f  ^  •  ♦  ***  rtttfttfl  1j7t*I  *r  Tl  !r*tf  n:  r  t < 

{ A', ^ f^4^+ttHtttH^rrt •  r * H *rHi-UrH 


i, 


V" 


•  •  «  *\4,  ,  ♦. 


•  •»* 

*  1  .'  i!  •  • 


.  - r^.x.  rn-n+rmH-}4iHt^tt  •  H+tffliU 

‘*‘w  — :  w^*'^..'*^."jt'T:ftH4'1'^i''ti^,'^i''‘ 

>,-•“■•  v4 /^  ■  ^  v, f' '•  ,^-h  t . '■'■'T 

•.*’.%,  A.  A  A  t^f  *"~*| 

'v‘  ij  V  V  ' 

.  .  ftCAlf 

•  '  *  j  i  . , ,  i  .  I  :  rit*  1  ’  *!’;.»  i  M  R »  ■ «  •  '  ?  ■  t  •  •  f  i  <  •  '  ’  t  W 
\  ~, 


DEPTH,  ft 


TEST  1*C 


VELOCITY,  fp* 


4000  8000  1 2000  1 6000  20000  24000 


Figure  20.  Velocity  versus  Depth,  Test  1C 


39 


1 


0 


100 


r 


120  c 


M0  r 


ItV  r 

r 

1*<0  p 

i 

r 

200  p 


220 


240  ^ 


4000 


TEST  2  *  C  * 


VELOCITY,  FPS 


K000 


1 2000  1 ^000  20000 
f -  r  t  -  i  —  I —  T 


O 


$> 


* 

$ 


t 

$ 


❖ 

<s> 


Figure  22.  Velocity  versus  Depth,  Test  2C 


41 


*  4 


COMMON  DEPTH  CROSS  HOLE  TEST 


y 


1*1  I  f  U  ‘  u  v 

ft 


»0O 


'  4  r  * 


148 


15  r 


1  54 


180 


103 


106 


189 


17  2 


1  75 


1  70' 


181 


104 


V 


4,  4  4  A  a,  a 

v  v  y  v  y 

■  j-  +  '  s> 

/  ly  v  *- 


V' 


w 


?• 

> 


T*V 


k  k 


,n  i-  M  A  A  ^  ^  *. 

,)  ’  ,  V  v;  ^  ^  •' 


m  A  sk  A  A  A 

'  y  ,»  y.  'J  v‘  v 


'A* 


>4*.  A  A  A.  * 

4  V  1  <  *  A  ”t 


&  A  4  A  A  v 

*•,  V  w  v"  • 


«» i  >  t  i  ,i  i  *  •/>  A 

v*.  v?  v  >;  v  y  w  vv*  - 


n,^.„  Jl  A  .  vftk.  . 1 4^  j -»  ,. 

i-‘  •*  v  v-  y,  ^ 


hr 


^  A  A  ^ 

\/  '✓•/  t ;  ^ 

A  A 


-*■" - V_r*>  .V 


A  >  A  •..  j«n,  .<►, 

^  v;  y  ''  v  "  " 

„  4  .4  A  .4  1  ^  A.  m.  ~  - 

■J  v  •;/  ,>  L:  •'  '-  •'•  •.••* v  v'V 


4 


•  *  - 


hi  V 

A  L.  A  A 


/V  V^v.  -  <•*», 


k*  V  'V  w 


••  !,!*.■*.  '.^.-—A,  /V.^, 

vf  v’ 


k  /  V  j| 


***' 


l  4,  +  *  **  *  u»J  14  440 

•  mm 


!  i  ■,  i  !  i  :■!••  ii,  mi.i  j  ■■  K < • .  .  - 1  ■ 1 : ,  .  !  > •  -.  t  .t 


but  t ho  ( unne  1  was  approx  im  it  e  I  v  half  wav  between  source  and  receiver  as 
opposed  to  the  othei  test  where  the  location  of  the  tunnel  is  closer  to 
either  the  source  or  the  receiver.  The  tabulation  of  arrival  time,  distance, 
and  associated  velocities  are  shown  in  Table  6,  and  a  plot  of  velocity  as 
a  function  of  depth  is  graphically  displayed  in  Figure  24.  The  time  domain 
signal  shows  not  only  time  delay  in  the  tunnel  zone  but  a  loss  in  amplitude 
and  h i gh  frequency  signal  content. 

41.  Test  ‘SC.  Observing  Figure  II.  it  will  be  seen  that  Test  5  was 
conducted  between  the  same  borings  as  Test  4.  The  test  configuration  was 
different,  however,  in  that  the  source  was  maintained  at  a  constant  depth 
of  162  feet  and  the  receiver  moved  in  increments  of  3  feet  in  borehole  C 

to  create  a  fun-tvpr  effect.  The  signal  quality  observed  in  Figure  25  (time 
domain  as  a  function  of  depth)  is  quite  good  from  a  depth  of  142  feet  to 
184  feet  and  the  tests  conducted  at  a  depth  of  163  feet,  which  is  common 
with  Test  4C,  exhibit  a  comparable  signature.  The  tabulation  of  arrival 
times  and  computed  velocities  are  shown  in  Table  7,  and  the  graphic  display 
of  these  velocities  as  a  function  of  depth  is  Figure  26.  Although  it  will 
bo  observed  that  the  general  pattern  of  velocities  shows  a  decrease  in  the 
area  where  the  tunnel  is  located,  the  fan-type  test  does  not  appear  to  be 
as  definitive  as  either  the  common  depth  or  the  offset  depth  tests  for  this 
case. 

S1E  records 

42.  Data  acquired  using  the  analog  SIE  seismograph  and  associated 
FM  tape  recorder  were  analyzed  in  greater  detail  than  the  previous  data 
acquired  with  the  EG&G  unit.  The  previous  section  on  data  reduction  discussed 
in  some  detail  the  procedures  used  to  present  the  data  in  a  manner  favorable 
for  in-depth  analysis.  This  discussion  of  in-depth  analysis  contains  only 
that  data  which  was  considered  represent  at i ve  of  the  entire  test  series 

at  the  Idaho  Springs  site.  Only  those  methods  which  were  thought  to  have 
given  a  positive  indication  of  the  tunnel's  presence  are  presented  herein. 

43.  Fourier  ( FFT )  analysis  of  the  seismic  wave  train  provided  a  great 
deal  of  additional  information  concerning  the  ident 1 f f cat  Ion  and  delineation 
of  tunnels  located  between  boreholes.  It  was  demonstrated  that  the  FFT 
analysis  provides  data  which  a  1 1 ow  the  determination  of  the  tunnel  location 
from  the  frequency  and  amplitude  effects  of  the  wave  train  due  to  a  tunnel- 
seismic  signal  interaction.  Different  portions  of  the  tunnel  cause  different 
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I/’.’,  1/"',  1/8,  anil  181  feet.  Those  four  travel  paths  have  the  longest  appar¬ 

ent  t  i  ivi  t  i  me  through  the  tunnel  .mil  should  demonstrate  t  he  greatest  seismic 
wave  tunnel  interaction  effects.  As  can  be  seen  in  Figure  10,  in  all  four 
i.ises  the  pi  etium  i  nan  t  signal  strength  ftom  1 1 0  to  TOO  Hz  is  greatly  attenu- 
1 1  ■  A  -  t  i  ,■ 1 1 1  .  i.oi.i  !  strength  buildup  i  an  he  seen  between  100  and  220  Hz  in 

■•hots  1  ’  >  .  1/8,  181,  and  184  feet.  All  of  these  last  four  travel  paths 
»ii  I  In  mu'll  t  fie  lover  half  of  the  tunnel  where  an  energy  buildup  was  docu¬ 
mented  in  this  treipiencv  range  for  Tests  IF  and  11F.  The  lack  of  a  larger 

'  ■ ,  i '  ’  1  e  se  '.’Wei  t  I  e-jiieni  i  es  mu'.  In-  due  l  e  the  inclination  of  seismic 

wave  I  l  .ivi  I  path  in  relationship  to  the  tunnel.  In  Test  SC  the  seismic 
wave  tiain  of  shots  1/S,  1/8,  181,  and  184  feet  actually  have  an  apparent 
travel  path  lor  a  significant  distance  in  areas  above  and  below  the  tunnel 
which  we t e  shown  in  IF  and  1  IF  to  highly  attenuate  the  seismic  signal. 

Hue  to  this  situation,  a  large  frequency  shift  n.iv  not  he*  achieved. 

4h.  Test  13F  (Frequency  Domain).  The  results  of  Test  1 3 F  are  similar 
to  Test  IF  (Figure  11)  except  for  the  shorter  travel  path  and  the  tunnel 
location  which  is  closer  to  the  transmitter  than  the  receiver.  In  shots 

144  ,  147,  ISO,  1 S  3 ,  159,  and  180  feet,  an  unaffected  signal  is  recorded 

wit*,  halt  ot  tiie  energy  buildup  at  frequencies  from  3  50  to  450  Hz  (Figure 
t.’l.  these  t  requeue  i  es  in  Test  13  F  are  higher  than  Test  3F  due  to  the  shorter 
travel  distance  in  3F.  Although  this  is  only  a  small  dispersive  effect 
over  a  short  distance,  it  will  greatly  increase  with  longer  travel  times. 

Such  phenomena  need  to  be  carefully  studied  before  an  undocumented  tunnel 
is  investigated  in  the  frequency  domain  at  larger  drill  hole  spacings. 

In  shots  at  140  and  15b  feet,  weak  coupling  between  the  drill  hole  wall  and 
tin-  seismic  source  or  the  geophone  results  in  a  poor  signal  strength.  Shots 
it  182,  1 n 5 ,  and  IbH  feet  demonstrate  the  same  signal  attenuation  above  the 

tunnel,  as  documented  in  3F.  As  before,  most  of  the  energy  above  220  Hz 
has  been  filtered  out.  At  171  and  174  feet  (ray  paths  which  transect  the 
lower  fra  1  f  of  the  tunnel)  a  large  energy  shift  from  higher  to  lower  frequen¬ 
cies  Is  evident.  This  is  similar  to  case  3F,  but  in  this  instance  a  much 
larger  portion  of  the  energy  is  transformed  from  higher  to  lower  frequen¬ 
cies.  This  frequency  shift  appears  to  be  more  than  100  Hz  with  almost  all 
of  the  energy  in  the  signal  now  being  under  270  Hz.  Test  13F  also  has  one 
travel  path  (177  feet)  which  passes  near  the  tunnel  floor  where  the  signal 


52 


Fourier  Analysis  and  Signal  Energy  Summation  Plots  for  Test  5  C 
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Fourier  Analysis  and  Signal  Energy  Summation  Plots  for  Test  13  C 
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is  .il  !  cim.il  f<l  01  I  1  1 1  i-ri'd  in  the  same  wav  js  above-  the  roof.  Tins  indicates 
•  >  '.imil.it,  Iml  not  as  extensile,  s  t  less  / 1  r.ic  t  ur  l  ng  tt-gime  in  the  tunnel 
1  lool  . 

.7.  IV;,  t  1  IF  t 1  i  me  Domain).  Fin  j>i  imai".  purpose  ot  this  if  i  sruss  i  on 
i '•  to  i  1  1  ust  t  ,tt  e  the  significance  ot  i  dent  1 1  i  ah  1  e  events  wti  i  c  h  influence 
the  entile  wave  signature  at  the  receiver  over  a  time  domain  of  approximately 
Mi  milliseconds  after  initiation  of  the  seismic  pulse.  In  addition  to  the 
(list  at  rival  of  the  I’-wave ,  which  under  ordinary  circumstances  will  be 
the  shot  test  t t ave  l  path,  othet  travel  paths  will  he  discussed.  The  presence 
ot  an  air-tilled  anomaly  (tunnel)  ictween  the  two  boreholes  creates  the 
possibility  tor  multiple  reflection,  refraction,  and  diffraction,  and  some 
ot  these  possibilities  will  be  addressed. 

id.  One  secondary  source  thought  to  be  a  primary  influencing  factor 
on  tlie  overall  wave  train  is  the  reflection  that  occurs  when  the  initial 
wave  encounters  the  air-tilled  tunnel.  The  length  ot  its  travel  path  will 
obviously  be  related  to  the  position  of  the  source  and  receiver  relative  to 
the  fixed  position  ot  l lu  tunnel.  In  the  conduct  of  Test  1  IF,  the  sequence 
began  with  the  source  and  receiver  at  a  common  depth  well  above  the  top 
of  the  tunnel  As  the  source  and  receiver  were  lowered  in  the  boreholes, 
the  travel  path  ot  the  reflected  wave  became  progr e ss i ve  1  y  shorter  until 
it  was  almost  equal  to  the  length  of  the  direct  wave  path.  This  occurred 
with  the  source  at  a  depth  of  lib  feet  and  the  receiver  at  a  depth  of  lb!> 
feet  (the  top  and  bottom  of  the  tunnel  are  at  depths  ol  approx imat e 1 y  161 
feet  and  17.'  feet,  respectively). 

49.  A  schematic  of  the  test  concept  is  shown  in  Figure  33.  Note 
that  in  addition  to  the  st ra i glit - 1 ine  travel  paths  between  source  and  receiver 
the  possible  paths  of  the  reflected  signal  from  the  top  of  the  tunnel  are 
also  shown  when  the  source  and  receiver  were  positioned  hove  the  tunnel. 

A.  third  event,  the  arrival  of  the  reflected  signal  from  the  ground  surface, 
will  also  he  considered.  Although  numerous  other  possibilities  tor  reflec¬ 
tions  proh.ihlv  exist,  the  discussion  will  tie  concluded  with  the*  diffracted 
signal  which  occurs  when  the  tunnel  lies  between  the  source  and  receiver. 

‘id.  A  tabulation  sutnmar  i  z  ing  depth  of  source  and  receiver,  distance 
to  the  receiver  by  way  of  the  direct  travel  path,  a  reflected  and  diffracted 
signal  from  the  tunnel,  and  a  reflected  signal  from  the  ground  surface  is 
compiled  in  Table  8.  The  corresponding  P-wave  travel  times  related  to  each 
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of  these  I'Vi'nt  s  .no  also  shown  in  the  tabulation.  This  tabulation  will 
bo  usfd  in  conjunction  with  the  time  domain  signal  to  discuss  various  events 
which  can  he  observed  on  the  seismic  signature. 

'il.  Figure  14  shows  the  digitized  time  domain  signal  recorded  by 
the  vertical  receiver  over  a  period  of  approximately  50  mi  1 1  i.scwonds  in 
each  of  the  14  shots.  The  first  event  which  can  be  readily  observed  is 
the  arrival  of  f he  first  signal  which  is  the  I’ -wave  train  journeying  along 
the  shortest  direct  travel  path  between  the  seismic  source  and  the  receiver. 

A  tabulation  of  these  travel  times  appears  in  column  7  of  Table  8.  When 
the  tunnel  begins  to  affect  the  length  of  the  travel  path  (Records  10-13), 
it  can  be  observed  that  the  signal  is  diffracted  and  delayed  when  compared 
to  the  uninterrupted  direct  travel  path.  A  maximum  time  of  8. 5  msec,  a 
delay  of  about  4  msec,  will  be  observed  on  Record  17.  In  this  case,  the 
locations  of  source  and  receiver  are  almost  on  an  Imaginary  centerline  through 
the  central  part  of  the  tunnel.  The  presence  of  the  tunnel  also  affects 
the  signature  of  Recoids  1-6  which  are  located  well  above  the  top  of  the 
tunnel.  A  high  degree  of  distortion  in  the  ela.-sed  time  frame  of  5  to  12 
msec  is  caused  by  t  fie  signal  reflected  from  the  top  of  the  tunnel  arriving 
about  1  to  2  msec  later  than  the  direct  wave.  Compare  columns  7  and  8  of 
Table  8.  This  distortion  can  be  used  as  a  possible  tunnel  identifier. 

52.  A  possible  diffracted  path  was  also  computed  using  an  average 
velocity  determined  from  the  direct  path  measurements.  As  an  example,  shot 
no.  14  was  recorded  with  the  source  at  a  depth  of  171  feet  and  the  receiver 
at  a  depth  of  180  feet.  The  direct  travel  path  was  47.5  feet  and  a  possible 
secondary  path  would  have  been  around  the  top  of  the  tunnel  resulting  in 
an  overall  distance  of  about  56.0  feet.  Its  travel  time  would  have  been 
about  1  millisecond  longer  than  that  of  the  direct  path.  However,  no  notice¬ 
able  effect  from  this  theoretical  possibility  is  observed  on  Record  14. 

This  diffracted  arrival  either  did  not  exist  or  its  signal  amplitude  was 
too  small  to  cause  a  visible  distortion. 

i  i .  Observing  Records  1  and  14,  it  can  be  seen  that  a  primary  reflec¬ 
tion  (circled)  occurs  at  a  time  of  27.5  msec  on  Record  1  and  at  a  time  of  34.2 
msec  on  Record  14.  Using  these  times  in  conjunction  with  appropriate  source- 
surf  ace- rece i ver  distances,  an  average  vertical  P-wave  velocity  of  10,250 
fps  can  be  computed.  Consequently,  these  reflections  could  very  well  be 
the  product  of  the  P-wave  signal  traveling  to  the  ground  surface,  reflecting 
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from  that  interlace  and  bat  k  down  to  t  lie  receiver.  Using  this  velocity, 
which  was  determined  1  tom  two  h  i  gh-qua  1  1 1  y  reflections,  and  the  known  total 
distance  Sot  Records  2-1),  the  predicted  arrival  time  for  each  was  computed 
and  tabulated  in  column  9.  In  most  instances,  the  stir f ace-ref lec ted  signal 
tan  be  recognized  on  the  records;  however,  in  some  instances,  the  arrival 
of  other  events  tends  to  mask  the  reflection. 

''A.  Another  interesting  observation  can  be  made  concerning  the  presence 
of  the  tunnel.  During  t  fie  construction  process,  in  which  drilling  and  blast¬ 
ing  v.-re  nsfil,  t  great  deal  ot  fracturing  occurred  in  the  rock  around  the 
tunnel  (possibly  as  much  as  one  tunnel  radius)  (a  feet).  This  observation  was 
confirmed  hv  the  lower  velocities  measured  in  this  region  and  supports  the 
observation  made  by  King,  et  .  al.,  (1984).  King  stated  that  the  construction 
effect  extended  more  than  h  feet  from  a  tunnel  driven  through  basaltic  rock. 
The  existence  of  these  fractures  will  also  cause  severe  attenuation  of  short- 
wavelength,  high-frequency  signals.  Consequently,  low  frequency  signals 
would  be  expected  to  dominate  the  seismic  signature,  as  discussed  in  more 
detail  in  the  previous  frequency  domain  analysis.  This  effect  can  be  readily 
observed  on  Records  11  and  12.  Additionally,  Record  12  shows  near-perfect 
in-phase  merging,  of  t  fie  diffracted  signals  which  traveled  around  the  top 
and  bottom  of  the  tunnel  re-uniting  on  the  far  side.  This  can  be  explained 
by  the  fact  that  the  travel  paths  are  virtually  equal  around  the  top  and 
bottom  of  the  tunnel.  Record  11,  however,  shows  distortion  caused  by- 
slight  lv  out-of-phase  signals  which  traveled  around  the  top  and  bottom  of 
the  tunnel.  In  t hi s  case,  the  travel  path  is  appreciably  longer  around 
the  bottom  of  the  tunnel  than  around  the  top  of  the  tunnel.  Consequently, 
the  time  domain  signal  readily  shows  the  distortion  caused  by  the  slightly 
later  wave  train  lagging  and  combining  with  the  first  signal. 

V>.  In  most  cases,  the  S-wave  arrival  was  very  difficult  to  determine. 
Of, serving  Records  1-4,  a  distinct  wave  train  arrival  can  be  seen  occurring 
in  the  12  to  13  msec  range.  This  wave  train,  even  though  somewhat  distorted 
by  the  first  primary  F’-wave  reflection  from  the  tunnel,  is  most  likely  the 
S-wave.  Its  velocity  is  about  3,800  fps.  As  testing  neared  the  tunnel, 
the  S-wave  arrival  became  incoherent.  Record  14,  however,  displays  a  promi¬ 
nent  wave  train  arrival  at  an  elapsed  time  of  18  msec.  This  signature  can 
be  traced  and  overlain  on  Record  1  to  establish  a  correlation.  Since  the 
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